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£SJ , Abstract 

We obtain the RPA summed one-meson exchange potential between nucleons in symmetric nuclear 
matter at zero temperature, from a model which includes p, a, u and n mesons. The behavior of rho mesons 
inside the medium is first discussed using different schemes to extract a finite contribution from the vacuum 
polarization. These schemes give qualitatively different results for the in-medium rho mass. The results are 
discussed in connection with the non-renormalizability of the model. We next study the modified potential 
as density increases. In the intermediate distance range, it is qualitatively modified by matter and vacuum 
, effects. In the long-distance range (r > 2 fm), one observes the presence of oscillations, which are not present 

in free-space. Features on this distance range are insensitive to the renormalization scheme. 

m ' 
o 

O ' 1 Introduction 

00 ' 

0^ ' Quantum Hadrodynamics (QHD) designates a class of models in which nuclear interactions are described 
through effective relativistic Lagrangians of nucleons coupled to different kinds of mesons. Such models have 
been successfully used during the past decades to study very different situations, such as nucleon-nucleon scat- 
tering processes in free space [|l] or the nuclear many-body problem || (for a recent review of QHD models, see 
||). fn vacuum, the one-boson exchange approximation gives a reasonable approximation to nucleon-nucleon 
scattering data. 

The situation becomes more complicated in dense nuclear matter since higher-order diagrams have to be 
taken into account. It is known for example that short range correlations arise due to the hard repulsive core in 
the free potential. An important achievement of the late eighties was the reproduction of saturation properties by 
parameter- free, (relativistic) Dirac-Brueckner G-matrix calculations [Q, ^J. In this approach, repeated exchange 
of free mesons give rise to an effective potential with a smoothened repulsive core. On the other hand, it is 
also well-known that effects coming from the Random Phase Approximation (RPA) will modify the two-particle 
interaction at high densities. 

For example, in QED the lowest-order interaction potential between two static charges is the usual Coulomb 
potential. Inside a plasma, RPA effects modify the photon propagator, and correspondingly one obtains a 
screened Debye potential. At larger distances, new phenomena appear. If the temperature is sufficiently low, 
the potential becomes oscillatory and damped as some power-law of the distance. These are the so-called 
Friedel oscillations and are originated by the sharp profile of the Fermi surface at low temperatures. Friedel 
oscillations are supported by a large body of experimental evidence in metallic alloys (see e.g. Q). 

Similar screening effects are encountered in a QCD quark-gluon plasma when the quark-antiquark potential 
is calculated within the same approximations || . Such effects appear also if one considers the spatial dependence 
of static meson correlation functions at finite baryon density in the Nambu - Jona - Lasinio model || . In a 
nuclear plasma, the presence of the medium will also give rise to a modified interaction. The modifications of 
the potential due to RPA have been investigated so far in the case of the one-pion exchange jlO| and for the 
Walecka a - to model [[ll], |l2|] . In both cases, a screened potential was obtained that differs from the one in 
vacuum . At long distances (r > 2 fm), the potential shows also an oscillatory behavior, as in the QED and 
QCD cases. 
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Our claim is that RPA may give important effects, both quantitative and qualitative, and has to be considered 
in many-body calculations of nuclear matter. As a matter of fact, it was shown that RPA corrections give 
rise to a modification of the in-medium cross section of the same order of magnitude as Brueckner ones, and 
compatible with experimental data available so far. Other examples are the reduction of neutrino opacities 
in neutron star matter from RPA corrections of the NN interaction |l4|, and the measurement of the 
electromagnetic response function in quasielastic electron scattering experiments Q. 

Of course, both kinds of effects should eventually be considered. Within a diagrammatic perturbation 
expansion, they appear as summations of two distinct subsets of diagrams, namely particle- hole loops for the 
RPA, and ladder summation for the G-matrix mentioned above. Both are essential ingredients, since RPA 
summation is needed in order to describe long-range correlations, and, on the other hand, ladder summations 
are needed to create short-range correlations. 

It would be desirable to develop some approach to incorporate both kind of summations in a systematic way. 
As discussed in Wjl , however, this can not be done simply by using the RPA modified potential as a driving 
term for ladder summation. Such an approach would lead to inconsistencies and double counting, and one has 
to look for more elaborated techniques. Parquet resummation JlTj is the minimal extension that can accomplish 
this purpose, but its extension to relativistic strongly interacting systems still encounters technical difficulties. 
Besides theoretical problems, parquet techniques applied to interactions with a simpler structure, (scalar X(j) or 
QED) already lead to very complex results exceeding standard computational resources and tricky convergence 
issues. [] 

Before embarking on such an attempt, our aim will be more modest. We think that it is useful to check 
whether the RPA effects give rise to sizeable corrections in the density and momentum transfer ranges of 
interest. Actually, we already know from other contexts that the answer is positive jl3], [H| ^J. Here we focus 
on the analysis of the modifications of the one-meson exchange potential due to RPA corrections to the meson 
propagator inside nuclear matter. We will pay special attention to the comparison of the modified potential, 
in contrast to the free-space potential. Keeping in mind the above considerations, since the RPA summation 
is not the only effect on the many-body problem, it cannot be used in naive ladder calculations. We insist, 
however, that this summation is a necessary (and important) piece of the problem, and deserves further study. 
In the course of the calculations, we obtain the modified meson propagator inside the medium, a result which 
is interesting by itself, since it allows us to study the effect of the medium on the meson propagation, the 
modification of the meson mass, the appearance of new branches, etc. 

In this paper, we concentrate on these modifications within a model which includes p mesons, in addition to 
the 7r, a and uj mesons already considered in the previous references. The inclusion of p mesons is a necessary 
ingredient in determining the proper isospin dependence of the nucleon-nucleon potential within meson exchange 
models Also, obtaining the rho propagator in dense matter will allow us to investigate how its in- medium 
mass changes as density increases. This issue has become important in the context of dilepton excess at low 
invariant masses observed in heavy ion collisions p3]. As we will show, however, the evolution of the p meson 
mass with density in this QHD model is very sensitive to the way in which the vacuum terms of the polarization 
are treated. This is due to the non-renormalizability of the derivative part of the p meson coupling which was 
introduced for phenomenological reasons. While renormalizability might not be a requirement in the case of 
effective models, one has to provide a prescription to eliminate the divergences that arise in these relativistic 
models consistently. As a matter of fact, one finds that, in renormalizablc QHD models, vacuum effects eliminate 
the pathologies which appear when all vacuum terms are simply left out (23). We will examine here several 
different prescriptions, leading to different consequences for the rho mass in nuclear matter. 

This paper is organized as follows. In section 2 we define the model and calculate the corresponding RPA 
meson propagators in the medium. In section 3 we discuss how to extract a finite vacuum term from the vacuum 
polarization in the case of rho mesons. We give in section 4 general expressions for the one-boson exchange 
potential obtained from the modified propagator. In section 5 we discuss the rho-meson dispersion relations 
and effective mass in connection with the renormalization procedures discussed in section 3. Results concerning 
the in-medium potential are showed in section 6. We first consider only the exchange of rho mesons. After 
this, we give some plots which are obtained by adding all kind of mesons present in our model. Mechanisms 
which could be responsible for smearing away the Friedel oscillations present in this potential are discussed 
in Section 7. Finally, our main results are summarized and commented in Section 8. We have gathered in 
Appendix A and B the explicit formulae concerning the matter and vacuum part of the polarization tensor of 
rho mesons respectively. Appendix C contains the configuration space contributions from different mesons to 

1 Interesting attempts in this direction may also be found in a non relativistic formalism in Jl8| [u], . 
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the nucleon-nucleon potential. 



2 Meson propagators in symmetric nuclear matter. 

In this section, we describe the main steps necessary to obtain the meson propagators in the RPA approximation. 
The starting point will be a QHD model. For this, we adopt the one described by the following Lagrangian : 

C = C N + L a + C u + C« + C p + Ci + C cc (1) 

which assumes nucleons interacting with several kinds of mesons : cr,n,p and to. In the latter equation, Cn 
corresponds to the nucleon Dirac free Lagrangian: 



where ip = ( ^ p I is the nucleon isospin-doublet (protons and neutrons) field with mass m, 



C N = (i/2) [ ifyy- d ip - (d tp) -7 r/>] - m tj) i> (2) 

1 i« t.Vip nnHpf 

C a = (l/2)[(d» *)(d v a) -p a 2 a 2 } (3) 

= (l/2)[(d» 7T)(d u 7?) -/^ 2 7f 2 ] (4) 

C p = -(1/2) [(1/2) ftp ■ ft mv - Mp 2 p ■ M (5) 
C u = -(1/2) [(1/2) FJT ■ F^ v - /xj 2 u " ■ w v ] (6) 

are the free Lagrangians for the mesons, with masses pi (i=cr, tt, p and u>). Ci gives the meson- nucleon interaction. 
We have adopted simple Yukawa couplings for the a, tt and uj mesons. For the p-meson we added a tensor term 
ft). Therefore, we have : 

Ci = g a 4> a tJj + 4> - i g-x 4> ~f 5 t ■ ir 

+ g p i>^r- - A ^ ^ ^ ? .Q v p - (7) 

with the following notations : 

= c^w 1 ' - 0" (8) 

i?/ 1 ' = p v - a* (9) 

Finally, £cc contains the counterterms necessary to eliminate the divergences from the vacuum polarization. 
We will discuss this topic in more detail in the next section . 

From Eq. (@) one can derive the equations of motion for the nucleons and mesons. In order to obtain the 
meson propagators, we have used the linear response theory around a given ground state of the meson-nucleon 
plasma. Our formalism is based on the introduction of Wigner functions. This formalism has been described 
in several papers in connection with the nucleon-nucleon interaction in a nuclear medium [ fi"l[ p3| , p4j| . We will 
give here only the main steps relevant to the model adopted above, and refer the reader to these papers (and 
references therein) for a detailed description of the method. 

The Wigner function for the nucleons is defined by : 

F(x,p) =< F[x,p) >=tr p E F(x,p) (10) 

where F(x,p) represents the nucleon Wigner operator 

H^P) = 7A4 / d 4 Re^ R ?(x + f ) ® V(*-f) (11) 

and pe is the equilibrium density matrix operator. For a system at a given temperature T allowing for particle 
number variation, this is given by the Grand Canonical operator. Given a quantum operator O, the statistical 
average is defined, as in Eq. ([To|), by 

< 6 >= tr \p E 6] (12) 
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(the symbol tr means the trace with respect to the quantum states available to the system). 

We assume that equilibrium can be described by the Hartree approximation. Within this approximation, 
meson fields are treated as classical fields and replaced by their statistical averages. These are furthermore 
restricted by symmetry properties. For symmetric nuclear matter, only < a > and the time-like component 
< lu° > survive [|j (all other meson mean fields are zero). This means that the nuclear background is described 
by the Walecka model. The mean-field values have to be obtained sclf-consistcntly, by solving the implicit 
equations 

Ma < a > = 9*Tr J d 4 p F H (p) (13) 

/£ < ^ > = -9u Tr J d 4 p j°F H (p) (14) 

Here, Tr stands for the spin-isospin trace. In a spin saturated system, parity conservation implies < tt >= 0. 
The value of the p mean field is related to the difference between neutron and proton densities, and since we 
are working in symmetric nuclear matter, we have < p >= 0. 

The nucleon Wigner function in the Hartree approximation is given by 

F H (P) = (7-P + M) (*f f °A (15) 

The nucleon effective mass is M — m — g a < a >, and P^ is defined by P' 1 = p^ + g u < > . The 
last matrix in Eq. ( |i~5| ) corresponds to the isospin structure of Fh{p) , and f(p) is the relativistic distribution 
function of nucleons, defined by 

fip) = (2^ 6{P2 ~ [Q+iPo) + n ~ iPo) ~ ff(_Po)] (16) 

where H(x) is the Heaviside step function and Q + (p ) (J7 _ (p )) are the nucleon (anti-nucleon) occupation 
numbers: 

n±G*) = h(±p ) [l + J (P0T , )} (17) 

In the latter equation, p is the nucleon chemical potential and /3 = l/T (we take the Boltzmann constant 
k B = l). 

The next step in the linear response formalism is made by introducing small perturbations of the meson 
fields and the nucleon Wigner function around their Hartree values. Therefore, one has to replace : 



<j(x) —* < a > +6a(x) 

w"(a;) -> < uj" > +Suj"(x) 

7?(x) — > 5 n(x) 

p»(x) 6p>*(x) 

F(x,p) - F H (p) + 5F(x,p) 



(18) 



in the corresponding equations of motion, and consider only terms which are linear in the perturbations (notice 
that < 7? >=< p ^ >= , as discussed above). One then obtains a system of coupled equations which can be 
solved for 6F(x,p) . When substituted in the remaining equations, it gives a system of equations for the meson 
sector. For symmetric nuclear matter, the resulting system is considerably simplified. First, the pions and the 
rho mesons decouple from the rest. Only the a and u> mesons are mixed together. However, the corresponding 
a + lu equations are the same that appear in the Walecka model ]Tl] , p5[ . Similarly, one reproduces for pions the 
well-known results ^ f23|. For this reason, we will concentrate here on the rho mesons. In this case, because of 
isospin symmetry, the equations of motion are the same for 

D^{k) b Pv = + ( k 2 



and p , and can be written as 

+ n^(fc) ] 8 Pv = o 



(19) 



Here, n'J I/ (fc) is the p-meson polarization tensor due to NN particle- hole loops. Explicitly, 



p 



d p Tr 



f, 



k, 



(7-(P -'o)+ M )(9 P r + ^ a k a )(l-(P ) + M)(g P Y ~ ^ vP k ) 



fiP 



2 

k/2) 



2m 



2m 



fip + k/2) 



p.k 



(20) 
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(Analytical formulae are given in Appendices A and B). This tensor can be decomposed into two terms (we 
omit Lorentz indices): 

u p (k) = n7'(fc) + n; ac (k) (21) 

The matter polarization, n™ a '(fc) vanishes in free space, i.e. at zero density and temperature. The second, 
vacuum term H v p ac (k) , gives a non-zero contribution even in free space. It contains divergent integrals from 
which one has to extract a finite contribution. This will be discussed in the next section. 

The p-meson propagator in the medium G{k) can be obtained by analyzing the response of the meson field 
against external nuclear sources. This procedure has been described in detail in . As showed in this reference, 
the resulting propagator is given, in a matrix form, by 

G(k) = -[D(k)]- 1 (22) 

and obeys the Dyson equation (also written matricially) . This equation, graphically represented in Fig. 1, reads 

G(k) = G°(k) + G"(k)U p (k)G(k) (23) 



where G°(k) is the non- interacting p-meson propagator in vacuum. Eq. (|l9|) can be analyzed with more detail 
by choosing a particular frame with the z-axis along the direction of k , in such a way that k^ — (w,0, 0, q). 
Then, the only non- vanishing components of the polarization are U° p °(k), H p A (k) = H 22 (k) ,H Q p 3 (k) — H 30 (k) and 
Ilfik). Eq. |[|) then becomes 



D 00 (k) 



D 30 (k) 





C n (fc) 








D 22 (k) 




D 03 (k) 



D 33 (k) 





Spo(k) 




8pi(k) 


X 


5 P 2(k) 




L 5 P3 (k) J 



(24) 



where 



D 00 (k) = 

D 11 (k) =D 22 (k) = 

D 33 (k) = 

D Q3 (k) =D 30 (k) = 



2 ) 



u 2 

( k 2 p 2 p n^(fc) ) 



„2 ,.2 TTll 



-w 2 + nl + nf{k) 



-wq 



U 03 (k) 



(25) 
(26) 
(27) 
(28) 



The dispersion relations can now be obtained by equating to zero the determinant of Eq. (|24| 
be factored out into two terms : 



which can 



\D{k)\ = (D"(k)) 



D 00 (k) D 03 (k) 
D 30 (k) D 33 (k) 







(29) 



We obtain two equations, corresponding to transverse and time-longitudinal modes. These are given, respec- 
tively, by : 



K( k ) 







(30) 



and 



k 2 p 2 - k 2 /q 2 n°°(fc) = (31) 
It is convenient to give these results in a fully covariant form. The p meson polarization can be decomposed 



as 



rfr\ v 
rf 



with 



k^k 1 - 



rf 



rf = 



—k» 

k 2k 



(32) 



where is the 4-velocity of the medium. H p t and H p l are Lorentz scalars representing the transverse and 
longitudinal components of the polarization. The dispersion relations read 



2 



n 



pT 



o 



n 



pL 



o 



(33) 
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The propagator is obtained by inversion of the dispersion relation, and is given in the general case by 

1 k^k v 1 1 



^ - ~GpL^ - G pT l , ^ 2 fc2 



GpL — 



GpT — 



In the referential where the fluid is at rest — (1, 0), we have the following relations: 

n pi = - k 2 / q 2 n° p °(k) ; n pT = - n"(fc) 



(34) 



(35) 



3 Vacuum polarization 

In this section we discuss the procedure to extract a finite contribution from the vacuum polarization tensor. 
In the case of the a and u>, or it mesons with a pseudoscalar coupling, the Lagrangian is renormalizable. This 
means that divergences can be eliminated consistently, at all orders, by introducing appropriate counterterms 
in the Lagrangian and imposing some physical conditions. This has been discussed in several papers; for a 
discussion within the formalism of Wigner functions see 23, |24|, ^6|. As pointed out in these references, the 



vacuum polarization gives an important contribution and eliminates some of the pathologies that appear in 
the semi- classical approximation (when all vacuum terms are simply discarded^]). Therefore, one would like to 
conserve vacuum effects. Here, however, there is an important difference when considering the rho mesons. As a 
matter of fact, due to the derivative coupling appearing in Eq. (0), the Lagrangian becomes non-renormalizable, 
which implies that the counterterms needed to compensate the infinities at a given order of approximation will 
not be valid at higher orders. Nevertheless, it is still possible to eliminate the divergences at a given order by 
the procedure described below. 

The method proceeds by dimensional regularization and the introduction of a counterterm Lagrangian. 

The vacuum contribution arises from the last term —H(—Pq) in Eq. (fi6|). Since the only tensors at our 
disposal in vacuum are g^ v and k^k v , and H p l " vac should be orthogonal to k^k v , it must be of the form: 



where 



npVen(^) 



Q^(k) = 



4tt 



k^k v 



9" 



After performing the trace, a straightforward calculation leads to 



n„ 



=(*) 



16 g 



it) 3 3(2tt) 3V ' 3\2mJ (2tt) 3 



3 (2 

8Mfc 4 ( /, 
W) 



2 fc 4 



3 \2mJ (2tt) 3 



2m 



(36) 
(37) 

(8M 2 + k 2 )I 2 
(38) 



The integrals X\ and X 2 diverge. These divergences can be extracted by the procedure of dimensional regular- 
ization 



Ti = 



al 4 pS(p 2 -M 2 )H(- Po ) 
d 4 p 



, S(p 2 -M 2 )H(- Po ) 



{p.k) 2 -fc 4 /4 



ll e9 = irM 2 



2 ~~fc2" 



1 1 f M K 

— + In — Ai 

e \ m 



_i + l n ^A 2 )+^(fc 2 ,M 2 ) 
e \ m 



(39) 
(40) 



9(k 2 ,M 2 ) is a known finite function given in Appendix B. Here, e is an infinitesimal quantity and Ai, A 2 are 
arbitrary finite constants. Therefore, there appear the following divergences in n p vac when e — > oo: 



1 



in the g p term , 



— in the -^—g p term, 
e 2m 



M 2 , k 2 . , 
and — m the 

e e 



fp_ 
2m 



term 



(41) 



2 Moreover, the vacuum contribution depends on the plasma thermodynamical state, so that it can not be simply subtracted in 
a fully consistent way. 
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At the one-loop level, it is possible to compensate the divergences by adding the following counterterms to the 
Lagrangian: 

C cc = {(A + B a + C o 2 ) RT ■ R MV + D (d a ■ (d a R^)} (42) 

in which A, B, C, D are constants to be determined by the renormalization procedure. In order to obtain the 
finite (renormalized) vacuum polarization, we have followed the same subtraction scheme as we used for the 
other mesons |23| . Namely we first obtain the rho field equation from the generalized Euler equation^ 

3C \ „ „ ( QC ^ „ 



d^ g ( 

dp v M \3(c^p„) 



\d(d a dnp u ) 



= 



with the full Lagrangian including the counterterms. After linearizing and Fourier transforming, the dispersion 
equation (|l^) reads 

[-k»k v + (k 2 - ifysT + K" + {A + Bo + Co 2 + Dk 2 ){k^k v - k 2 g^)] 5p v = (44) 

Now the vacuum rho polarization can be made finite if A, B, C, D consist of an infinite part cancelling the 
divergences we have pointed out in Eq. (41), and a finite part which we now proceed to determine on the basis 
of physical arguments. We will have 

(A + Ba + Ca 2 + Dk 2 ){rk h 

2M 2 



p ren 



p vac 



9p 
3tt 2 



m 



k 2 



1 



2 



AY J_ 

2m / 6ir 2 



QM 2 \n— + 8M 2 9 + k 2 (lii — 
m m 



2M / /. 



2m 



with the decomposition 

„2 



A 
B 



3tt 2 



-- + lnA 2 
e 















f e 




m 






1" 


2m 


( 9 P fp 


~ 2 




\ 2m 



(a + (3 a + 7 a 2 + 5k 2 ) j {fc V - ^k"} 



(45) 



h In Kn 



fj _^9a_f 9pfp 

1 ir 2 \ 2m 



D 



A x 


\ 2 9l 


3 


2m, 


I 3tt 2 


e 


A' 


i'-L 


1 


2m y 


' 6tt 2 


e 



1 , A 1 

h In A 2 

e 2 



A 

2m 



~^ 2 ~ 



P \ 2 9 



2m / 37r 2 
4 lnA 2 - InAi - 2 



- + 4 In A 2 

e 



InAi 



(46) 



In order to determine the finite constants a, /?, 7, (5, some physical constraints have to be imposed. The standard 
p5[ requirements are that: 



* (1) The vacuum polarization has to vanish in free space, when the rho meson is on its mass shell (here, 
"shell" ={k 2 = p 2 pl M = m,o- = 0}): 



n 



p vac\ she 



n = 0, 



(47) 



* (2) The field equation can be written in the standard form by defining effective coupling constants <? c ff and 
/eg. Requiring that the effective coupling be equal to the bare one in vacuum leads to the condition 

u±1 o vac 



dk 2 I shell 







(48) 



* (3) The new couplings aR^R^ , a 2 R^ V R^ V , that we had to introduce in order to subtract the infinities, 
should not appear in the dispersion relation on the rho meson mass shell, leading to the conditions 



L p vac q 

da (shell 



(49) 



3 The D counterterm makes this generalization necessary 
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and 
*( 4 ) 

dtt 



p vac 



(50) 

This amounts to four relations to determine the four constants a, (3, 7, 8. 

Alternative sets of conditions can be imposed, leading to different expressions of the vacuum polarization. 
In particular, we can recover the renormalization scheme suggested by Shiomi and Hatsuda ]27j as a special case 
of our general formalism. We have also tried these alternative schemes in our calculations in order to investigate 



the influence of the treatment of the vacuum. As we will see in section (5.3), the various schemes predict a 
completely different behavior of the in-medium rho meson mass as density increases. 

The peculiarities of each scheme, outline of derivation and corresponding analytical expressions are given 
in Appendix B. They can be subdivided into two classes: The first class of schemes comprises those which are 
designed not to introduce new couplings with the a mesons (conditions (3) and (4) of the method outlined above). 
This is the case of the scheme presented here, labelled "scheme 1" and of a very similar one labelled "scheme 2" . 
With a different point of renormalization, the scheme presented in the case of the a meson in the classical paper 
of Kurasawa and Suzuki @] also belongs to this class. We applied it to the p meson and labelled the result 
"scheme 6" . The second class of renormalization schemes preserves the structure of the regularized expression, 
introducing counterterms a + /3er + 7c 2 so that they can be factorized into aM 2 /m 2 = a(m — g a a) 2 /m 2 . As a 
result, conditions (3) and (4) cannot be fulfilled anymore. We called the scheme obtained in this way "scheme 
3". The method of Shiomi and Hatsuda (scheme 5) and a related one used by Sarkar et al. (29) ("scheme 4") 
belong to this class. 

The advantages and caveats of the respective renormalization schemes will be further discussed in a forth- 
coming workQ. 

4 One-boson exchange potential in the medium. 

In this section we calculate the one-boson exchange potential obtained after RPA summation. As mentioned 
in the previous section, when the background is symmetric nuclear matter, rho mesons (and pions as well) 
decouple from the other mesons, whereas a and to mesons are coupled together. Therefore, the total potential 
shows the following structure : 

V = V a+ ^ + V* + V 71 (51) 

where V a+U) ,V* and V™ are obtained by a + u>, one-pion and one-rho exchange, respectively. Our method 
follows a similar procedure to the construction of the one-boson exchange potential in vacuum |Q . The essential 
difference is that free meson propagators are replaced by in- medium propagators. Also, the external nucleon 
lines correspond to in-medium spinors, as they arise by solving the nucleon Dirac equation in the Hartree 
approximation. It is important to mention, however, that one arrives to the above result by computing the 
energy associated to a pair of interacting nucleons inside the plasma, a result which ensures that no double 
counting has been made. The method has been given in detail in pT| , and applied to the calculation of V a + M . 

As in the free case, we expand the resulting potential in powers of jj, where p represents the momenta of 

the external nucleons, and keep only terms of the order ~ {~mY ■ This has the advantage to simplify the spin 
structure of the nucleon-nucleon potential. For degenerate nuclear matter, one has p ~ pp, where pp is the 
nucleon Fermi momentum. If the density is close to saturation, the next term in the expansion would be a small 
correction of the order ~ (fpp) • Of course, if density is much larger, higher order terms in the above 

expansion should be taken into account. 

In constructing the one-rho potential one has to take into account the diagram shown in Fig. 2, where 
the double line represents the in-medium rho-meson propagator, given by Eqs. (|2^ - ^8|) and Eq. (|34|). The 
amplitude corresponding to this diagram is given by 

M p = - 



SabG^(k) |[x 2 ^(p' 2 ,4)] (i9pl u ~^" a W ~P2)^j r\ Mp 2 , S2 ) X2 ]} (52) 
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In this equation, Xi > X2 stand for the quantum states of the initial nucleons with four-momenta pi and p2- They 
have associated spinors u(pi,si) and u(p2,S2) ■ Similarly, the prime indicates the corresponding magnitudes 
for final nucleons. The isospin operators rf , t\ (a, & = +,—, 0) label the associated rho charged fields. From Eq. 
([52]) we obtain the one-rho exchange potential, by using the procedure outlined above. In the center-of-mass 
frame (CM) of the nucleons^, the potential is given by: 



VMQ) = 



-G pL 1 



+ 4M<? P 



+ AM 



1>L 
2m , 



2M 2 8M 2 2M 2 



(?A Q).S 



2 1 fp 

2m 



-GpL 

G p t 



AM 2 2M 2 
i 



AM 2 



2M 2 



— (qAQ).S+(a 1 .a 2 ) 



AM 2 



(a 1 .q){a2-q) 
AM 2 



2 r 



G p T (cfi.fa) 



q 2 {ai.q){a 2 -q) 



AM 2 



AM 2 



T\ ■ T2 



(53) 



In the latter equation, <7i (fj) represent the spin (isospin) Pauli matrices for the two (i — 1,2) interacting 
nucleons. In the chosen frame, initial nucleons have momenta p and — p , whereas final nucleons are assumed 
to have momenta p ' and -p ' . We have introduced the notations : 



q=p-p 
G pL (0,q)- 



(p + p')/2 



•ft 



n°°(o,9) 



G pT (0,<t) 



q 2 + f l 2 p + m/(0,q) 



(54) 



Next we will construct the in-medium (RPA summed) nuclear potential in configuration space. This is done 
by Fourier transformation of the corresponding momentum-space magnitudes. Before proceeding, we include a 
phenomenological form factor in the nucleon-meson vertices. This amounts to making the replacements : 



9a 

F a {q) 



g a ■ F a {q) 
fp ■ F P (<i) 

+ q 2 



(55) 



(a = a,u>,Tt,p) in all previous equations in momentum space. After some algebra, the configuration space 
potential can be decomposed in the following way |l| : 



V(f) = V c (r) - - {V 2 V NL (f) + V NL (f)V 2 ) 

+ V LS (r) L-S + V ss {r) #i ■ <r 2 + V T (r) S 12 



(56) 



with the following notations : L is the angular momentum operator, S = h (<Ji + $2) is the total spin of the 
nucleons and the tensor operator. In this way, V c [r) is the central part of the potential, Vls(t) is the 
spin-orbit, Vss(t) is the spin-spin, and Vr(r) is the tensor potential. In addition to the spin structure, one has 
to include a factor f\ ■ T2 in the pion and rho contributions (omitted here). The second term in Eq. ( |56| ) is 
non-local. Its contribution will depend on the explicit form of the two-nucleon wave function and, therefore, 
can not be directly plotted. For our analysis in the next section we will not consider this term. 

Each one of the pieces in Eq. ( |56"| ) has contributions from different meson exchanges (a + uj sector, 7r or p). 
Their expressions can be found in Appendix C. 

4 We consider here the case where the rest frame of the background fluid coincides with the center of mass of the collision. For 
the effects due to a relative velocity between both frames, the reader is referred to ]3l] 
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5 Results for the rho-meson propagation 



5.1 Choice of parameters 

We will analyze the behavior of the in-medium potential as density changes, at zero temperature. Our calcu- 
lations include the vacuum and matter polarization in the meson propagators. Since the vacuum contribution 
does not vanish at finite qin general, we cannot take directly the parameters available in the litterature. Instead, 
we must perform a fit anew for each one of renormalization schemes described in Appendix B. We have fitted 
the value of the coupling constants g a , Qui, 9p, fp and cutoff parameters A CT , A w , A w , A p in such a way that, at 
zero density, the model reproduces as close as possible the Bonn 'Potential B' [Q (table A. 3 in this reference). 
The value of the pion coupling was kept constant since it is well known from experimental data. The values 
of meson masses are taken directly from this reference. However, we do not include the 5 and rj mesons in the 
fitting procedure. Since the Bonn potential has been adjusted to reproduce low-energy nucleon scattering data, 
we expect that our model can give a reasonable description to these experimental data. We have chosen to 
make this fit in the region from 0.5 to 2.5 fm for all mesons. The resulting parameters are given in Table 1, and 
the values of the obtained chi-squared for the different pieces of the potential appear in Table 2. 





9<T 


9w 


9p 


fp/dp 


A(j 


A w 


A. 


Ap 


Bonn B 


10.6 


17.55 


3.36 


6.1 


1900 


1850 


1700 


1850 


set 1A 


8.10 


16.267 


4.216 


4.974 


1630.55 


1562.19 


1179.48 


1881.93 


set IB 


7.55 


15.817 


3.390 


5.115 


1558.22 


1550.12 


729.40 


1365.54 


set 2A 


7.83 


15.783 


3.718 


4.996 


1390.92 


1493.97 


976.30 


1797.30 


set 2B 


7.59 


16.091 


3.882 


4.319 


1532.25 


1521.51 


690.48 


1112.22 


set 3A 


8.54 


15.910 


3.862 


5.811 


1080.65 


1441.00 


1212.30 


1638.73 


set 3B 


8.53 


16.377 


5.710 


4.547 


1329.54 


1500.82 


1150.89 


1260.05 


set 4A 


7.22 


14.624 


3.303 


7.033 


1540.38 


1629.70 


1144.72 


1375.11 


set 6A 


7.57 


14.809 


3.636 


5.731 


1626.54 


1590.76 


1125.17 


1711.54 



Table 1: Values the coupling constants and cutoffs which adjust the Bonn potential B |l| for each renormalization 
scheme 





X\V C ) 


X 2 (Vss) 


X 2 (Vt) 


X 2 (Vls) 


set 1A 


0.096 


0.026 


0.014 


0.078 


set IB 


0.136 


0.232 


0.019 


0.179 


set 2A 


0.078 


0.049 


0.056 


0.172 


set 2B 


0.123 


0.182 


0.061 


0.203 


set 3A 


0.026 


0.039 


0.035 


0.142 


set 3B 


0.038 


0.065 


0.031 


0.082 


set 4A 


0.079 


0.079 


0.054 


0.210 


set 6A 


0.087 


0.034 


0.069 


0.191 



Table 2: \ 2 values corresponding to the parameter sets given in Table 1, for each component of the NN 
potential: central (C), spin-spin (SS), tensor (T) and spin-orbit (LS) 

In these tables, the number labelling each set is determined by the corresponding renormalization scheme 
(see Appendix B). The letter A means that the best fit for all components of the potential was looked for. As 
will be seen in the next section, the dispersion relation of the p meson may display heavy meson modes and, 
more annoyingly, zero-sound modes for too high values of the p meson cutoff. Therefore we performed the 
fit once again, but now with the constraint that a lower cutoff is used for the p meson in order to avoid the 
appearance of these spurious branches. The letter B corresponds such fits, and due to this restriction they come 
with slightly higher \ 2 values. 
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We do not give parameter sets for renormalization scheme 5, because this scheme was constructed by Shiomi 
and Hatsuda in such a way that the vacuum polarization vanishes identically (for all q) in the vacuum. In this 
case, the expression of Machleidt Q is recovered exactly, therefore there is no need of readjusting the parameters 
and the original values of Machleidt (Bonn potential B) may be used. 



5.2 Rho meson dispersion relations 

With the above values of the meson parameters, we have first performed a numerical study of the rho-mesorj^] 
dispersion relations defined by Eqs. (|30|), (31). Vacuum appears through the renormalization scheme discussed 
in Section 3. We use the renormalization scheme described in that section (scheme 1 of Appendix B), as well 
as alternative schemes described in Appendix B. A sample of the resulting branches are plotted in Figures 3 
to 5 for longitudinal and transverse modes, for various values of the nucleon Fermi momentum pp = 0.3, 0.4, 
0.5 (all magnitudes are given in units of the nucleon bare mass m). In our model, this correspond to densities 
Po, 2.37 po and 4.63 pa with po being the nuclear saturation density. The variables u> and q are defined by 
k = (u, 0, 0, q), as in Section 2. 

We at once observe that the dispersion relation strongly depends on the chosen renormalization scheme, 
both directly, and indirectly since it conditions the choice of parameters at the time of fitting the potential 
(see Table 1). Let us first discuss the general features of the transverse mode for parameter sets 1A to 6A, 
which were obtained by fitting the NN potential without imposing restrictions on the allowed range for each 
parameter. 

In the timelike region, besides the normal branch, we also have two heavy meson branches in the general 
case. This kind of meson branches have been obtained in other models of QHD, and originate from vacuum 
effects ]23| , p4[ . Such meson modes are present at densities corresponding to pp > 0.3 for all renormalization 
schemes, with the exception of parameter set 4A, corresponding to the renormalization scheme of Sarkar ]29| ]. 
It must be noted, however, that this parameter set has a value of the p meson cutoff appreciably lower than 
other A-sets. As will be discussed later, a smaller value of the cutoff has the effect of reducing or eliminating 
the heavy meson branches. 

At high density, there appears moreover a zero sound branch in the spacelike region. This latter branch 
could have especially strong (and unpleasant) effects on our results for the potential, since it appears as a pole 
in the propagator. We found such a branch for renormalization schemes belonging to both classes considered in 
this work: in the "increasing rho mass" class, with parameter set 2 A at pp > 0.4 and also for the "decreasing 
rho mass" renormalization scheme of Shiomi and Hatsuda at pp > 0.6. Nevertheless, we will see in the sequel 
that this pole can be eliminated by reducing the value of the p meson cutoff. 

The longitudinal modes have a similar structure in the timelike region, with a normal branch and two heavy 
meson branches. In the spacelike region no zero sound branches were found. 

In order to eliminate the spurious branches, it was possible to find a second series of parameter sets with 
stronger cutoffs for the p meson, at the cost of a slightly degraded quality of the fit. We called the resulting 
sets IB, 2B and 3B. This permits to eliminate the zero sound branch at all densities investigated (up to pp=2) 
and to reduce or eliminate the heavy meson branches. It was moreover checked that the choice of a lower cutoff 
parameter does not affect appreciably the position of the normal branch. 

The dispersion relations for parameter set IB are shown in Figure 3. A moderate cutoff A p =1365 MeV 
was applied. The zero sound branch is now removed, but we still have heavy meson branches. The left panel 
displays the transverse modes. The right panel compares the transverse (thin line) to the longitudinal (thick 
line) modes for the normal branch. At q = the transverse and longitudinal modes coincide; for finite q they 
differ only slightly. Finally, we note that the intercept of the normal branch with the q = axis goes to higher 
frequencies with increasing pp (or densities). This is related to the fact that the effective mass of the p meson 
increases with this renormalization schemes, as discussed in next section. 

The dispersion relations for parameter set 2B are shown in Figure 4, now with a strong cutoff A p =1112 MeV. 
Only the normal branch remains for the transverse mode. The longitudinal modes are also appreciably cleaned, 
with only a persistent heavy meson mode at high momentum transfer. A still lower cutoff would eliminate it, 
but would spoil the accuracy of the fit of the NN potential. 

In Figure 5 we show the transverse dispersion relation obtained with a renormalization scheme of the second 
("decreasing rho mass") class. The parameter set used is 3B, with a reasonable p cutoff A p =1260 MeV, resulting 
in a clean dispersion relation: there only remains the normal branch. Moreover, the quality of the fit remains 

5 The analysis of the a , ui and n dispersion relations has been performed in [ pj| |. 
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very good. We do not show the longitudinal modes since we found that longitudinal and transverse modes are 
indistinguishable. We note that the intercept of the normal branch with the q = axis is lower at finite density 
than in vacuum, and slightly increases from pp=0.3 to pf=0.5, corresponding to the fact that the rho meson 
effective mass with this renormalization scheme first decreases with density and reaches a minimum around 
Pf = 0.35 (see next section). 

We also programmed and plotted the dispersion relations obtained with schemes 4 to 6, but chose not to 
display them here . Scheme 6 ("Kurasawa-Suzuki") was found to be essentially similar to other schemes of the 
same class (1 and 2) but more problematic when we tried to remove the sound mode and heavy meson branches. 
For the dispersion relation obtained with the scheme 5 of Shiomi and Hatsuda |^7j with the parameters of 
Machleidt Bonn potential B, we have again a normal branch and two heavy meson branches in the longitudinal 
as well as in the transverse modes. Moreover, a tiny zero sound branch also appears in the transverse mode 
at high enough densities (p > 7 po, pf <~ 0.6). We performed a fit with a lower value of A p and found a 
parameter set which behaves essentially like case 3B. For the scheme of Sarkar el al. (scheme 4) we obtain only 
the normal branch for both transverse and longitudinal modes. This renormalization scheme is free from heavy 
meson branches and zero sound modes. This is due in part to the fact that the fit to the potential in vacuum 
requires in this case typically lower rho meson cutoffs of the order of 1300 MeV. Still we found this scheme 
more robust than the previous ones, since no such branches were found even for extreme values of parameters 
adjusting the potential. 

5.3 Rho meson effective mass 

From the dispersion relations one can obtain the in- medium mass p, e f / , defined as the solution p, e f f = u> of the 
dispersion relations Eq. (^), ([H]) at q = 0. The transverse and longitudinal modes yield the same effective 
mass, since the polarizations are equal in this limit. In Figs. 6-8 we show how p e ff evolves, in units of the 
vacuum mass /i p , as pf grows (remember that pp/m = 0.3 corresponds to saturation density). In particular, we 
wish to investigate the role of the vacuum on its behavior, in connection with the claim of Shiomi and Hatsuda 
[^7| that the vacuum contribution solves this issue. 

We first show in Figure 6 the effective mass obtained by keeping only the matter contribution in the rho 
meson polarization (thus neglecting the vacuum term). This figure was obtained with the parameter set of 
Machleidt potential B (see table 1). If we would limit ourselves to a moderate density and frequency, it would 
appear that the rho meson mass is slightly increasing. A closer inspection reveals, however, that the structure 
of the dispersion relation is more complex and has 3 solutions for pp < 0.246 (that is, a density p < O.bpo), 
corresponding to anomalous branches intermingled with the normal one. At higher density, only the higher 
branch survives. This type of behavior is in fact well documented. It was first observed in the case of the u> 
meson by Lim and Horowitz p2| and also appears in the case of the pion ]33]| . Contrarily to a rather widespread 
belief, it is not cured by the introduction of a stronger form factor. As a matter of fact, we redraw on the same 
figure the rho mass for a low cutoff (1 GeV) and see that the problem, although less acute, is still present. 
Moreover, we cannot choose the cutoff as we please, since it was obtained from a fit of the nucleon-nucleon 
potential in the vacuum. 

In the case of the a, uj and tt mesons, the situation is largely improved by introducing the vacuum term. 
There, a lower normal branch cleanly separates from higher heavy meson branches. Without any cutoff, two 
heavy meson branches appear at zero density and merge and disappear at high density. The heavy meson 
branches can thereafter be removed by the application of a reasonable cutoff. As we will see, the same occurs 
in the case of the rho meson. 

In figure 7, the on-shell renormalization schemes 1 and 2 of Appendix B were used. The parameters were 
chosen so as to obtain the best fit of the nucleon-nucleon potential in each case (parameter sets 1A and 2A 
of table 1). We also plotted on the same figure the resulted obtained by adapting the scheme of Kurasawa 
and Suzuki for the a, uj mesons to this case (parameter set 6A). In this figure and the next one, only the 
mass corresponding to the lower normal branch is represented. There also exist two branches at higher masses 
(typically 3.5 p p and 7 fi p ) which merge and disappear at very high density, which we do not show here. It is a 
common feature of all these renormalization schemes that the lower branch gives an in-medium mass which is 
larger than in vacuum. In all cases, the increase is quantitatively similar for all the three schemes, and somewhat 
stronger for scheme 6 of Kurasawa and Suzuki. On the other hand, the heavy meson branches largely depend 
on the value of the cutoff, so that this cutoff can be used to remove these branches altogether. As an example, 
we give in Table 1 a parameter set which makes it possible both to obtain a satisfactory fit of the potential and 
to remove the heavy meson branches (set 2B). 
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On the other hand, if one makes use of the renormalization scheme proposed by Shiomi and Hatsuda in 
[ p7| , one then obtains a rho mass which decreases with density. We have re-derived the formulae given in this 
reference using our own methods (scheme 5 of Appendix B) and checked that our result is in agreement with 
theirs. The resulting effective mass is plotted in Figure 8. As advertised, the rho meson mass decreases until 
Pf = 0.35 (or p = 1.6p ), and then slowly increases. We note that there also appear heavy meson branches. 
We plotted on the same figure the result obtained from a similar scheme used by Sarkar et al. |2^] (scheme 4 of 
Appendix B). While Shiomi and Hatsuda subtract the vacuum at all momenta k, Sarkar et al. set the vacuum 
to zero at k 2 = [i 2 in a way similar to our condition Eq. (f47|). The normal branch obtained from scheme 4 is 
quantitatively very similar to the result of Shiomi and Hatsuda. 

These results are tantalizing since they are in agreement with the currently accepted concept of a decreasing 
rho meson mass based on other models and experimental data. Now, in order to interpret the contradiction 
which appears between schemes 1,2,6 on one hand, and 4,5 on the other hand, we have developped a further 
scheme (scheme 3). As further explained in Appendix B, the introduction of counterterms (A + B a + C a ) 
does not preserve in general the original structure of the diverging term proportional to M 2 = (m — g a o} 2 . It 
can be required that this M 2 structure be preserved by dropping conditions (49,5^). The resulting effective 
mass is also plotted on Figure 8, and we see that it actually decreases as well. 

The idea of a decreasing rho mass in the medium was popularized by the Brown and Rho scaling conjecture 
|54| as a result of chiral symmetry restoration. It allows to parametrize the in-medium mass of vector mesons 
as 

^■wl- (0.18 ±0.05)— (57) 
Mp Po 

as a function of the density p (written in units of the nuclear saturation density po). For example, as a 
consequence of QCD sum rules, Hatsuda and Lee J35| or Leinweber and Jin J3(| obtain a decreasing mass. 
The use of QCD sum rules, however, was criticized for containing some uncertainties and inconsistencies with 
chiral perturbation theory |37|. Another caveat is the parametrization of the spectral function used in these 
calculations |38|, J39], which further reduces the predictive power of QCD sum rules arguments. Finally, recent 
QCD sum rules calculations with an improved vacuum subtraction p0[ further contributed mitigating the simple 
Brown-Rho decreasing mass picture. 

Other calculations, based on chiral approaches (see e.g. GJ, and j4^] for a review) or quark-meson coupling 
models |43j also predict a slight lowering of the rho meson mass. In chiral models including a pmr interac- 



tion piece (see also the brief discussion in this work, section |7.l| ), the rho mass is only slightly modified, the 



important effect there is the broadening of its spectral function. It is interesting to note that, depending on 
the renormalization scheme, the position of the maximum of the spectral function is shifted up- or downwards 
depending on the choice of the renormalization procedure (Q] vs. [^H)). The theoretical situation is far to be 
settled, and there is a possibility for the rho mass to increase, rather to decrease, due to medium effects (see 
also e.g. @, 0). 

There are a number of experiments in which it has been claimed that one can extract information about the 
vector mesons in a nuclear medium. In heavy-ion collisions experiments as HELIOS-3 |l£| and CERES Jig] , 
the excess of dilepton production at invariant masses lower than the bare rho meson mass might be explained 
in the framework of the vector dominance model by assuming a dropping of the in medium rho-meson mass 
[p2| [50|[5l[]. This interpretation is sometimes called, for short, the B/R (Brown-Rho) scenario. However, it has 
been pointed out that keeping a constant rho mass and introducing a medium modification of its width |5^, |5j| 
could also explain the dilepton excess. Friman and Pirner [Q have also argued that the contribution of higher 
resonances (*N(1720)) to the rho-meson self energy in matter plays an important role in shifting the strength 
to lower invariant masses. The 'broadening scenario' sometimes receives the name of R/W (Rapp-Wambach 
H|) scenario. 

For some time, measurements of polarization-transfer experiments with polarized protons seemed to favor 



the dropping-mass hypothesis as well 55 |56j, |57| ]. New analyses of the results with a better treatment of 
relativistic effects, however, concluded that the data could be explained without this assumption p^ |. 

In summary, the problem of how medium effects will modify the rho-meson properties is still not well settled, 
both from the theoretical and experimental points of view, although the possibility that its mass will be smaller 
at higher densities seems to be more favored. Future experiments, like HADES at GSI, will hopefully help 
clarifying this issue. 

In our hadronic model, we have seen that the inclusion of the vacuum contribution largely affects the meson 
propagation in matter. Vacuum renormalization still allows for a increasing, as well as for a decreasing meson 
mass, and does not by itself settle the issue. If we believe that the rho meson mass should be decreasing, then a 
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renormalization scheme which preserves the original dependence in the nucleon effective mass should be chosen. 
But then, in order to be coherent, the same procedure has to be used for the <r, 7r and lj mesons as well §0) 0. 

In view of this situation, we have made some comparisons using the two classes of renormalization schemes 
described above, each one of them giving a different behavior of the rho meson inside the nuclear medium. 



6 Results for the RPA effective potential in the medium. 

We shall now analyze the different components of the in-medium potential in position space, obtained in section 
4, as density changes at T = . As discussed in that section, the mixing terms of the p-meson polarization with 
other mesons vanish if one considers symmetric nuclear matter, a circumstance which will allow us to consider 
first the contribution to the potential of this meson alone. Later, we shall add the contributions of the a, lu and 
7r mesons. 

In performing the present analysis, it is useful to make the following division in the distance variable r, as 
was made in |ll| : 

a) In the short-range region < r < 1 fm, relativistic effects (beyond the quadratic approximation considered 
above) become essential. Also, in this distance range the extended structure of the nucleon must be properly 
taken into account. The introduction of form factors allows to simulate this effect for not-too-short distances, 
such as 0.5 < r < 0.8 fm. For shorter distances, calculations from our model will not be reliable. 

b) The intermediate region will be defined as the distance range 1 < r < 2 fm. Within this range, the 
vacuum potential becomes quantitatively modified by medium and vacuum effects, as compared to the free- 
space potential. 

c) In the long-range region : r > 2 fm, matter polarization dominates, and one can observe qualitatively 
new features. 



6.1 Rho-meson exchange 

On the left upper panel of Fig. 9, we show the contribution of one-rho meson exchange to the central component 
of the potential. We compare the results obtained for the in- medium potential at pp/in = 0.4 (corresponding to 
a density p — 2.37po) to the free-space potential (dotted curve), which was calculated setting pp — 0. Vacuum 
polarization effects are kept in all cases. We made the calculation for three types of renormalization schemes, one 
of the first "increasing rho mass" class (scheme 1 with parameter set IB), and two of the "decreasing rho mass" 
class (scheme 5=Shiomi-Hatsuda with Machleidt's Bonn B parameters and scheme 3=ours with parameter set 
3A). 

From this comparison, we observe that the use of a screened interaction in matter can bring appreciable 
modifications to the potential described by the exchange of free mesons in the vacuum. In particular, in the long- 
range region, the interaction is qualitatively different from vacuum. Instead of the usual exponential damping, 
at non-zero density the potential becomes oscillatory due to the presence of Friedel and Yukawa oscillations. 
This behavior is obtained whatever renormalization procedure is chosen. 

Friedel oscillations arise because the analytically-continued matter polarization shows branch cuts starting 
at q = i2pp in the complex g-plane. This phenomenon was first discovered for a QED plasma [|| and was 
also evidenced in the case of a QCD plasma Q. The corresponding analysis for a nucleon plasma in different 
models was performed in |uj [l^, [Oj. On the other hand, Yukawa oscillations appear when the analytically- 
continued boson propagator has a pole (in the complex q-plane) away from the real and imaginary axis. Such 
a phenomenon has been found so far for a nucleon plasma in the one-pion exchange approximation, and for 
a quark-gluon plasma when one-gluon exchange is considered p9| . In order to separate Friedel from Yukawa 
effects for rho mesons, one needs to perform the analytical continuation of the meson propagator and to study 
the evolution of the Yukawa pole as density and temperature evolves. This will be the subject of a future work. 

These oscillating phenomena can have consequences if one goes beyond the Hartree approximation, by 
including the polarization contribution into the ground-state energy. In the case of Friedel oscillations, it 
has been found that a periodic-density configuration, with period equal to the characteristic period of Friedel 
oscillations, has a lower energy than a constant-density configuration |B3|. This can be interpreted as a transition 
to a spatially-structured configuration. 

For a fixed renormalization scheme, we obtained that the repulsive core extends farther when density in- 
creases. On the other hand, for a given value of the density, the first minimum occurs earlier for renormalization 

6 Preliminary studies indicate however that this would spoil the behavior of the a and n effective masses [BOI] . 
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schemes of the first "increasing rho mass" class than for schemes of the second class. In the long-range zone, 
the period of oscillation is seen to be the same for all renormalization schemes, and depends only on the density, 
as expected from the physical origin of Friedel oscillations. 

Matter polarization dominates the potential qualitative features at large distances. The choice of a renor- 
malization scheme can, however, introduce some quantitative changes on these features. The amplitude of the 
oscillations depends indirectly on the renormalization scheme. For example, to each scheme, there corresponds 
a parameter set with different values of the couplings and cutoffs. The amplitude can be enhanced by the 
presence of a dip in the dispersion relation of the transverse modes, indicating the proximity of a zero sound 
mode, even when no pole actually appears. On the other hand we found that the appearance of such a branch 
depends on the renormalization procedure, the second class of renormalization schemes being less liable to the 
appearance of zero sound modes. Nevertheless, we would like to stress that the oscillations themselves are not 
a consequence of possible zero sound modes, since we found them also in the renormalization scheme of Sarkar, 
where no zero sound was obtained even for extreme choices of the parameters and density. 

Similar features can be observed in the remaining components (spin-spin, tensor and spin-orbit) of the 
potential. They are plotted in the three remaining panels of Fig. 9). They all show an oscillatory behavior. In 
all cases, the amplitude of oscillations increases with density. 

6.2 Combined-meson potential 

We will now present some results which are obtained by adding the contribution of the a + u> sector and the tt 
exchange to the potential given above^. For the central component there is no pion contribution, and one has : 

V c (r) = V° + "(r)+V£(r) (58) 

The result is plotted in Fig. 10 in vacuum (solid curve) and at finite densities pf = 0.3, 0.4, 0.5 (dashed, dot- 
dashed and dotted lines respectively). In this figure the renormalization scheme 1 was chosen with parameter 
set B. The free-space potential (which includes vacuum polarization with on-shell renormalization, scheme 1) 
has a potential well with a minimum at r ~ 1.5 fm. At finite density, the position of the first minimum is 
displaced towards shorter distances: at saturation density (pf = 0.3) it corresponds to r ~ 1.25 fm, whereas 
for pf = 0.5 it is located at r ~ 0.95 fm. There also appear secondary minima which depth can be significant. 
In the example chosen, at pp = 0.5 (p = 4.63 p$), the second minimum is situated at r ~ 2.1/m with a depth 
of -8 MeV. The right panel of Fig. 10 focusses on the long-range oscillatory behavior. Here it can be seen that 
the period of the oscillation decreases and its amplitude increases with increasing density. 

We also studied the spin-spin, tensor and spin-orbit components of the potential. The results are not shown 
here; basically the same features as already mentioned in the case of the rho component were observed. 

7 Mechanisms reducing the amplitude of oscillations 

In this section we discuss the mechanisms which could suppress the oscillations found in the in-medium NN 
potential. 

We will first investigate the contribution of meson loops to the p-meson self energy. As is well known, the 
spectral function of the rho meson acquires an important contribution of the pnn loop (see eg. [[44| ) . A meson 
with a mass distribution, such as the a or the p, could smear away some features in the potential by merging 
contributions with different ranges. As we will see, the effect on the potential is negligible. For consequences 
on the dilepton production in the vector dominance model, we refer the reader to the vast litterature devoted 
to this subject (see eg. [n], |62j and references therein). 

A second and more effective mechanism is the rounding off of the edge of the Fermi momentum distribution 
function introduced by short-range correlations or by a finite temperature. This will be investigated in the 
second part of this section. 

7.1 Meson loops 

We have up to now studied the contribution of nucleon-hole loops to the p meson polarization. On the other 
hand, the p meson also couples to the pion, a fact which is at the origin of the decay width of the p meson. Let 

7 The contribution from 8 and r\ mesons is not included here since they introduce only a small modification to the vacuum 
potential of [y 
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us add to the Lagrangian a piece 



W = gpir-K^TT X 7?).p M + ~g p ^{^ X 7T).(p M X 7?) 



(59) 



At the first order of the linear expansion method of pi] , meson loops do not appear. It would be necessary 
to go to higher orders in the cluster expansion, so that meson-meson correlators appear. We would like here 
to estimate the correction brought by the pnir polarization loop. From Green function techniques, one obtains 

M M M 



n 



pTTTT 



d 4 p 

(2tt)< 



(2p - kf{2p - k) v D v (p)D v (p - k) - 2ig- 



2 



d 4 P 



(60) 



In general, the pion propagator also includes medium effects. In order to keep things simple however, we will 
keep here the free pion propagator D®(k) — l/(k 2 — p 2 + ie). 

The coupling constant g p7r ir is determined by fixing the imaginary part of the polarization on the p meson 
mass shell and in the rest frame to its experimental value 



1 



= lmU p7r7T (u; = fi p ,q = 0) = 

p p 48irp p 



0*2 ~ 4,4)3/2 



(61) 



151.5 MeV, one finds g 2 /(^n) ~ 2.9 The real part of the polarization can be expressed in 



Putting T p7r ^ 
terms of the divergent integrals (M 



Tie IP" 



2 gt 



3 (2tt) 2 

After performing dimensional regularization. 



2T 1 +fc 2 (/4 



k^k 1 



(62) 



3(2^) 2 



g" 



k^k 1 - 



{—-£)9(tf,&)- — 



(63) 



The infinity is cancelled by a counterterm of the form A P T T7T R I1U R^, that is to say, the renormalization proceeds 
exactly as in the unambiguous case of the NN loop contribution to the u> meson. The finite result is: 



TZeW - 2 9/ " T7T 
3(2^) 2 



k^k v 



§0£ - k 2 ) (1 - 9(p%pD) + (p 2 It) (9(p%pl) 6(k 2 ,pD) 



LM; 



The imaginary part is given by 



n 2 / 2 \ 3/2 

p ™ 487r y 1 p 2 



9(k 2 -Apl) 



(64) 



(65) 



In the calculation of the potential, the imaginary part will not contribute since the polarizations entering 
the formulae are to be taken at to — while the imaginary part should be an odd function of u> so that the 
Onsager relations be fulfilled. The real part of the polarization yields a small shift to the dispersion relation. 
We have checked that its effect on the potential is negligible. At the level of approximation considered here, 
the shift is ^-dependent but does not depend on thermodynamical conditions. It would if we introduced the 
modification of the pion propagator by NA loops or a bath of thermal pions (see eg. |6^]). 



7.2 Effect of short range correlations 

A second mechanism which can attenuate the long range oscillatory behavior of the in-medium potential is the 
rounding-off of the momentum distribution by short-range correlations of the Brueckner type. As commented in 
the introduction, a fully consistent calculation at the level of parquet approximation would be required. Here, in 
order to estimate the order of magnitude of this effect, we will simply introduce such a momentum distribution 
instead of the Fermi-Dirac one in the calculation of the polarizations. 

Among the results available in the literature we chose two: a non relativistic Brueckner-Hartree-Fock cal- 
culation of Baldo et al. starting from a separable Paris interaction (6^], and a Dirac-Brueckner-Hartree-Fock 
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calculation by de Jong and Malfliet |65| . It seems more appropriate to use the calculation of de Jong and 
Malfliet since it is fully relativistic and was obtained starting from the Bonn potential. Moreover, these authors 
checked explicitly that the thermodynamical consistency be optimally fulfilled. On the other hand, comparison 
with experimental data seems to favor the larger depletion and smaller discontinuity found in nonrelativistic 
calculations such as that of Baldo et at, so that we also present results with this distribution. 

We found that the data published by de Jong and Malfliet was well reproduced by the following fit at 
saturation density: 

0.91- 0.055 ( — ) - 0.07( 1 )l°gf 1 | if k < p F 

\PfJ \ PF J \ PF J 



DBHF 



(k) = { 



0.035 



k/pi 



0.94 



exp 



k 

1.2 — 

PF 



0.3693 exp - 2.4 — 

PF 



if pf < k < 1.5pF 



if k > l.bpF 



(66) 



For the data of Baldo et al., we used 



' 0.79-0.13^ — V 0.19(l. -— Wfl-— 1 if k< PF 

\PF J \ PF J \ PF J 



F, 



NRBHF 



(k) 



16. exp - 4.5 — 

PF 



0.2127 exp - 1.6 



Pf 



if pf < k < 1.5pf 



if k > 1.5pp 



(67) 



We now calculate the polarizations at vanishing temperature as 



n 



BHF 



(k) 



tf 3 } 



an 



(T=0) 
FD 



(M 



dp 



F, 



BHF 



(p) 



(68) 



where H F d i s polarization which would have been obtained with a Fermi Dirac step function 0(pf — p) 
and is given in Appendix A. Next we introduce these polarizations in the expression of the potentials given 
in Appendix C. Figure 11 compares the central potential as obtained with parameter set IB at vanishing 
temperature and saturation density for the Fermi-Dirac momentum distribution (dashed line), the relativistic 
DBHF distribution of de Jong and Malfliet (full line) and the nonrelativistic BHF distribution of Baldo et al. 
(dotted line). It is seen that the oscillations are damped by Brueckner correlations as was to be expected. The 
effect is however not very severe and the oscillations persist. 



7.3 Effect of finite temperature 

A similar damping of the oscillations is produced by a non- vanishing temperature. We found that the Brueckner 
modification of the momentum distribution had the same effect as a temperature of ~ 10 MeV with a Fermi- 
Dirac distribution. It can be seen on Figure 12 that Friedel oscillations are rather robust against temperature 
and are not washed out for temperatures as high as 30 MeV. 

As a by-product, we obtained the behavior of the p effective mass as a function of temperature for all 



renormalization schemes. At saturation density, it was found that /U e // (calculated as in section ^3) slightly 
decreases with temperature for renormalization schemes of the first class 1, 2 and 6, whereas it slightly increases 
for renormalization schemes of the second class 3, 4 and 5. At four times the saturation density, p e ff was found 
to decrease for all renormalization schemes. The results are well described by a linear law for temperatures 
T G [0 - 100] MeV: 

^£l {p ,T) = ^(p,0)(l-a-) (69) 

with 
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a = 0.4 for p = pg renorm. scheme= 1 or 2 

a = 0.35 " " " " 6 

a = -0.12 " " " " 3, 4 or 5 



a = 1.5 for p = 4pQ renorm. scheme= 1 or 2 

a = 0.8 " " " " 6 

a = 0.55 " " " " 3, 4 or 5 



8 Discussion and Conclusions 

In this paper we have discussed RPA effects on the nucleon-nucleon potential obtained within the one-boson 
exchange approximation in symmetric nuclear matter at zero temperature. The model we used for our calcu- 
lations includes a , oj , ir and p mesons interacting with nucleons via Yukawa couplings. The Wigner function 
technique and linear response analysis were used to obtain the meson propagation in medium. The lowest order 
of this scheme is the mean-field approximation. When vacuum effects are properly renormalized, it becomes 
the Hartree approximation. The next order of this approximation has been shown to be equivalent to Green's 
function calculations of the meson propagators at one- loop order f23f. 

In our calculations, medium effects appear in two ways. First, in the nucleon legs of the one-boson exchange 
diagram the nucleon effective mass appears. Secondly, the meson propagators are calculated including nucleon- 
hole loops and vacuum polarization effects. 

It is legitimate to ask whether one should consider other loops involving mesons or resonances [Q. Only 
a brief discussion concerning pion-pion loops was given, since several thorough studies exist in the litterature 
p| ^33], [n| . Loops involving resonances can be expected to have a smaller but sizeable contribution of similar 
characteristics to the nucleon loop taken into account in this work. A serious study with renormalization would, 
however, go beyond the scope of this paper. 

In symmetric nuclear matter, rho mesons decouple from other mesons. Therefore, it is possible to study 
one-p exchange alone and simply add the contribution of er , ui and tt to the nucleon-nucleon potential at the 
end. We have first studied the dispersion relations which describe the propagation of rho-mesons in matter as 
density changes. One observes the presence of two types of branches : a normal branch, which is the analogous 
of the free-space mass-shell condition, and several heavy-meson branches. The latter appear as a consequence 
of vacuum effects, in all the QHD meson models that have been investigated |2^, |24|, |6(|. Here, due to the 
derivative coupling of the rho-mesons to the nucleons, the Lagrangian is non-renormalizable. Yet, it is possible 
to extract a finite contribution from the vacuum at each order in the cluster expansion^. This procedure, 
however, contains some arbitrariness in the case of non-renormalizable Lagrangians. 

In view of this, we have compared the results arising from two different classes of renormalization schemes. 
The first one is analogous to the one usually employed in previous works p5[ ^8, 10 and is designed to remove, 



as far as possible, the spurious new couplings to the sigma field which have to be introduced in the counterterm 
Lagrangian. The renormalization procedure suggested in p7f has also been used for comparison. We found that 
it belongs to a second class of procedures, where the new couplings to the sigma field are not minimized, but 
instead the original structure of the expression of the vacuum contribution is preserved. 

By using the first scheme, we find that the effective mass of rho-mesons grows with density, a feature which 
is disfavored by most of present theoretical approaches and by the interpretation of dilepton production data 
in heavy-ion collisions. The second method gives the opposite behavior : the in-medium rho-meson mass drops 
as density grows, more in agreement with present ideas. However, it is recognized that large uncertainties 
are contained both in the theoretical calculations |37j] and the analysis of experiments mentioned above 
f32[ [42] , p2[ . We have kept for these reasons these two schemes in our calculations. 

Next, we investigated the RPA nucleon-nucleon potential obtained by the exchange of one rho meson. At 
short and intermediate distances, it shows a repulsive behavior. In the long-range it contains new qualitative 
features, as compared to the free-space potential. It becomes oscillatory, with an amplitude which decreases 
with distance, due to the combination of Yukawa-like and Friedel oscillations. The appearance of this oscilla- 
tory behavior might give rise to a new phase of dense matter, characterized by a spatially structured density 
distribution. We have verified that the qualitative features within this distance range are not sensitive to the 
renormalization scheme used. One then expects that they are to be found in other models describing rho-meson 
propagation in dense nuclear matter. 

8 We use the word renormalization in a wide sense to describe this procedure, in spite of the non-renormalizability of the model. 
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When other mesons are added into the potential, we find a similar behavior. The potential is repulsive 
at short distances and attractive at intermediate distances. In the long-range, r > 2 fm, we find again an 
oscillatory behavior in all components of the potential, with the same origin as above. The amplitude of these 
oscillations increase with density. 

The robustness of these oscillations against several possible damping mechanisms has been tested. We first 
discussed the influence of pion-pion loops. In the spacelike zone explored by the potential, we found that it 
should not have a decisive impact on the results presented. This does not enter in contradiction with the large 
body of publications for which the spectral function is needed in the timelike region. 

We also investigated the amount of smoothing of the Friedel oscillations to be expected from modifications 
to the Fermi surface. In particular, we discussed the effect of the temperature and short-range correlations. 
We performed a calculation with a momentum distribution rounded off by short range correlations of the 
Brucckner type, and found they have the same consequences as a temperature of about 10 MeV, namely, that 
the oscillations experiment a moderate damping, and are seen discernible. They only disappear for temperatures 
higher than 30 MeV. 

One has to exert some care when considering the applicability to the above results to actual nuclear-matter 
calculations. The RPA summation concerns only a class of many-body diagrams, so that the addition of 
different diagrams, such as ladder diagrams, might modify our results (one would expect, however, that in the 
long-distance range the RPA results would hold). The obtained potential can not be used naively to perform 
more elaborated calculations (such as ladder resummation) since it would result in double counting of diagrams. 
Our results rather have to be interpreted as an indication of a new qualitative effect that one might find. 

To summarize, screening effects can appreciably modify the nature of nuclear interactions in a medium. 
Vacuum effects also compete to modify the interaction, so that it is important to design a procedure which 
incorporates them in a consistent way. However, the appearance of singular behaviors, as Friedel and Yukawa- 
like oscillations seems to be a rather general property of the nuclear interaction. 
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Appendix A: Polarization, matter part 



We give here the complete formulae for the matter part of the rho-meson polarization tensor in symmetric 
nuclear matter. The vacuum contribution is given in Appendix B. This tensor can be written as : 



n^(fc) = - J d 4 pi^( P ,k) j( P ,k) 



where 



I^(p,k) 



9p 



+ 



2m 



2m 



2 r 



9 P 



(4 M 2 + k 2 - 4p 2 ) g» v - 2 k^k v + %p^p v 

(4 M 2 k 2 + fc 4 + 4 p 2 k 2 - 8 (p.k) 2 ) g^ - (fc 2 + 4 M 2 - 4 p 2 ) k^k u 
+8 (p.k) (k^p v + p^k v ) - 8 fc W 



In this equation, M is the nucleon effective mass, and m its free mass. We have introduced : 

f(p + k/2)-f(p-k/2) 



J(p, fc) 



fc • p 



(1) 



(2) 



(3) 



The nucleon distribution functions, as they arise from the relativistic Hartree approximation, are the same for 
protons and neutrons (in symmetric nuclear matter). They are given by 

5{p 2 -M 2 ) \ H{p°) H(p°) 



[2nf 



e P(p a -»aff) + 1 + e l3(p"+P C ff) _|_ 1 



H{-p° 



(4) 
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Here, fi e ff = /i + <? w < w° > is the nuclcon effective chemical potential, /i is the true chemical potential and (3 
is the inverse temperature. The time-like mean field value < u>° > of the ui meson is also calculated within the 
Hartree approximation. Finally, H(x) is the Heaviside step function. 

At zero temperature, it is possible to perform the integrals in Eq. ([!]) analytically. Moreover, in order 
to obtain the nucleon-nucleon potential, within the approximations discussed in this work, we only need to 
calculate the matter polarization on the k° = axis. In this case, one gets more compact expressions. We give 
here these formulae. By defining q — \k\, e — y/ M 2 + (q/2) 2 and the nucleon Fermi energy Ep — \JM 2 + p 2 F , 
we can introduce the notations 



A{q) 



q - 2p F 

q + 2p F 



B(q) 



qE F + 2p F e 
qE F — 2p F s 



(5) 



One finds, after some algebra: 
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for the transverse part. 



Appendix B: Polarization, vacuum contribution 

The vacuum contribution arises from the —H(—po) term in the distribution function f(p). This contribution 
is divergent and has to be renormalized. After performing a dimensional regularization and subtracting the 
appropriate counterterms, we obtained in section III a finite expression for the vacuum polarization: 
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9a 



.95 



with the function 



9(k 2 ,M 2 ) = 9{y) = y 



dx 



i 



(i) 

(2) 



where y = 1 - k 2 /(AM 2 ). 

The choice of conditions to be imposed on the polarization to determine the finite constants a, (3, 7, 5 
determine various renormalization schemes. 

In order to shorten the notations, we define 
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J pram 
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(3) 



All these derivatives are to be evaluated in the point k 2 — /j, 2 , M — m. 

We will note 8 = 8(k 2 , M 2 ). Finally, we will work in units of the nucleon free mass, i.e. we set m = 1 in all 
expressions given in this Appendix. 

Scheme 1 

When the full set of conditions ([l7| - 15(]) is applied to the total regularized vacuum polarization of the p 
meson, the following expression is obtained 
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Scheme 2 



The former scheme has the drawback that we do not recover the expression of the vacuum polarization of 
the u) meson when we take the limit f p — > and replace g p by g u . The reason is that the structure of the 
infinities is not the same for the different couplings. By renormalizing all contributions together we are in fact 
introducing spurious counterterms where they need not be. For example, in the (f p /2m) 2 contribution we have 
both M 2 /e and k 2 /e, needing all sets of counterterms to cancel them, whereas in the g 2 contribution we have a 
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1/e which requires only the A R} w ■ Rp V counterterm. This can easily be cured by splitting the counterterms as 
follows A = A gg + Afg + Aff, B = Bf g + Bff, C = Cff, D = Dff, (where Af g , Bf g mean terms proportional 
to g p f p , and so on) and renormalized separately each contribution. We obtain 
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Scheme 3 

We observe that there appear infinities which can be written as depending on M (rather than on a), as for 
example M 2 /e. After absorbing the infinities, the finite part of the counterterms is ajj + bffcr + Cffa 2 . The 



constants a//, 6//, c// which are calculated by imposing the conditions ( |47| ) and ( |4S|j50| ) cannot be cast any 
more in a form proportional to M 2 . One would like to require that the M 2 structure be preserved by dropping 
conditions pgJSQ). Proceeding in this way, we obtain 
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Scheme 4 

We show in this section and the next how the results of Shiomi and Hatsuda [^?J and Sarkar et al. 
be recovered. We first observe that the regularized vacuum polarization can be written as follows: 
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If we^now impose that this structure be kept after the renormalization, with A and B substituted by finite 
A and B after extracting the infinities, we can determine the constants by imposing that A and B vanish on 
the mass shell. By so doing we obtain 
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and 
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We now use the fact that our function 9 is related to the integrals appearing in the expressions of Shiomi 
and Hatsuda and Sarkar et al.. 
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With these relations, it is straightforward to check that Eq. (g) coincides with the expression of Sarkar et al. 
Scheme 5 

The expression of Shiomi and Hatsuda can be recovered by the following recipe: First, one should subtract 
from the expressions of the functions A, B from the preceding paragraph their values in the true vacuum M = to, 
and then, replace them in the expression 9 Eq. ([?]). In terms of our 9 functions, one has 
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Scheme 6 

Another standard renormalization scheme is that used by Kurasawa and Suzuki p8| in the case of the a 
and u> mesons. In this scheme, the conditions Eqs. ([l7 48) are applied at the physical mass k 2 = fj, 2 whereas 
conditions (E^JsO) are taken at k 2 — 0. The following expression is obtained 
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9 We note that this is not simply the same as subtracting the true vacuum from the total polarization, since the effective mass 
appears outside of A, B in Eq. ([ij) 
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Appendix C: Interaction potential in spatial coordinates 

In this appendix, we give explicit expressions for the different meson contributions to the several pieces of the 
potential, as they appear in Eq. ( p6[ ) : central, spin-orbit, spin-spin, tensor, and the nonlocal potential Vnl(t)- 
A subscript will label these components, while the superscript indicates the meson which gives this contri- 
bution. 

These expressions are valid when the rest frame of the background fluid coincides with the center of mass 
of the interacting particles. 

Rho meson 



vy\ r ) = 



v$<r) 



V T "\r) 



V { L "kr) 



^ftoMrtt G PL (l- + 9 P (A) G pL 



~ 9 ? ° pT (4^) 



9 P + 2M 

9 P + 2M 

9 P + 2M 
1 

27r 2 M 2 r 3 
1 

2n 2 M 2 r 2 



-2 ji,r2 , 



2m 

A. 

2m 

A. 

2m 



1 f°° 
7r 2 M 2 - ^ dqq 3 sm(qr)F p G pT 

1 f°° 

^2 _/ o ^<Z 2 cos((zr)^ 2 G pT 



G 



pT 



8tt 2 M 2 



g sin(gr) f p i g p 



G 



pL 



dq q cos(gr) T p \ g p 
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GpL 



2^M l r jo 
with the following notations 



dq q sin{qr) T 



G 



pL 
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+ GpT 
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(1) 

(2) 
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(4) 
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-1 



^+^-11^(0,9) 



G pT = G p \q) = 



-1 



q 2 + V 2 p + n p \0,q) 



Sigma and omega mesons 



The mixed a-u> sector gives 



V™ (r) = J"" sin ^ G ° ( X + 4^2) - Guj L (l 

( f \ 



4M 2 



+9u Ft Gu) T 



I 4M 2 J 
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(6) 
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In the above formulae, one has 



dqq 2 cos(qr) T 2 G^t 



(7) 
(8) 
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G ff = 

G(JLJ 

Gto L = 

GuiT = 



g 2 +,i 2 -II° o (0,g) 



(g 2 + m 2 + n CT (0, g)) (g 2 + m 2 - II°°(0, g)) + (n° ff (0, g))^ 
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g 2 + M 2 +nii(0,g) 



Pi meson 



Finally, the one-pion exchange only contributes to the spin-spin and tensor components : 



where : 
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fig 
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(11) 
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Figure 3: Rho-meson dispersion relations for parameter set IB. The solid line correspond to vacuum density, 
and the dashed, dot-dashed and dotted lines to finite density with a value of Fermi momentum pp/m—0.3, 0.4, 
0.5 respectively. The left panel shows transverse modes. There are normal and heavy meson branches. The 
right panel compares longitudinal (thick lines) to transverse modes (thin lines) for the normal branch. We use 
the same conventions as in the left panel to represent various densities. Magnitudes are given in units of the 
free nucleon mass m. 
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Figure 4: Rho-meson dispersion relations for parameter set 2B. The left panel shows the transverse branch at 
various density with the same conventions as in Fig. 3. Only the normal branch remains in this case. The 
longitudinal branches are represented in the right panel. Besides the normal branch, we have a remnant of the 
heavy meson modes at high density. 
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Figure 5: Rho- meson dispersion relation for parameter set 3B. Only the transverse branch is shown. The 
longitudinal branch almost coincides with the transverse one so that it could not be distinguished by eye from 
the former. Only normal branches are present in this case 
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Figure 6: In- medium p- meson mass as a function of the nuclear Fermi momentum, when the vacuum term is 
discarded. For small densities, it would appear that the effective rho mass increases; however neglecting of 
vacuum fluctuations is at the origin of an unpleasant cusp, which cannot be removed even by very low values 
of the cutoff parameter in the form factor 
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Figure 7: In-medium p-meson mass as a function of the nuclear Fermi momentum. Vacuum polarization is 
included according to the first class of renormalization schemes 1, 2 and 6 (= "Kurasawa-Suzuki"). 
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Figure 8: In-medium p- meson mass as a function of the nuclear Fermi momentum. Vacuum polarization 
is included according to the second class of renormalization schemes 3, 4 {— " Sarkar") and 5 (= "Shiomi- 
Hatsuda" ) . 
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Figure 9: One-rho exchange component of the potential at 2.4 times saturation density (pp—OA) for various 
renormalization schemes: scheme 1 with parameter set IB (solid line), scheme 3 with parameter set 3 A (dashed 
line) and scheme 5 (= "Shiomi-Hatsuda" ) with the parameter set of Machleidt Bonn-B potential (dot-dashed 
line). The potential in vacuum is also shown for reference (dotted line). The left upper panel shows the central 
component, the right upper panel displays the spin-spin component, the left lower panel corresponds to the 
tensor part. Finally the spin-orbit component is represented in the right lower panel. 
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Figure 10: Central component of the potential including all (a, u, p, ir) mesons in the medium, using parameter 
set IB with renormalization scheme 1. We use the same conventions as in the left panel to represent various 
densities. The left panel represents the short range part of the potential. The right panel illustrates the 
oscillatory behavior in the long range part. 
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Figure 11: Influence of the rounding off of the momentum distribution function on the amplitude of Friedel 
oscillations. We show the central component of the potential including all mesons at saturation density and 
vanishing temperature, using parameter set IB with renormalization scheme 1. The solid line was obtained for 
a Fermi-Dirac distribution (as in previous figures). The dashed line uses the parametrization of the result of a 
(relativistic) Dirac-Brueckner-Hartree Fock calculation by de Jong and Malfiiet [58] . The dash-dotted line uses 
the parametrization of the result of a (nonrelativistic) Brueckner-Hartree Fock calculation by Baldo et al. [57] . 
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Figure 12: Influence of temperature on the amplitude of oscillations. We show the central component of the 
potential including all mesons at saturation density, using parameter set IB. The solid, dashed, dot-dashed and 
dotted lines correspond to temperatures T=0, 10, 30, 50 MeV respectively. 
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